Noise-enhanced trapping in chaotic scattering 
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We show that noise enhances the trapping of trajectories in scattering systems. In fully chaotic 
systems, the decay rate can decrease with increasing noise due to a generic mismatch between 
the noiseless escape rate and the value predicted by the Liouville measure of the exit set. In 
Hamiltonian systems with mixed phase space we show that noise leads to a slower algebraic decay 
due to trajectories performing a random walk inside Kolmogorov-Arnold-Moser islands. We argue 
that these noise-enhanced trapping mechanisms exist in most scattering systems and are likely to 
be dominant for small noise intensities, which is confirmed through a detailed investigation in the 
Henon map. Our results can be tested in fluid experiments, affect the fractal Weyl's law of quantum 
systems, and modify the estimations of chemical reaction rates based on phase-space transition state 
theory. 
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The responses of nonlinear dynamical systems to small 
uncorrelated random perturbations (noise) can be sur- 
prising and seemingly contradictory. Noise usually de- 
stroys fine structures of deterministic dynamics, e.g. it 
fattens fractals but it can also combine construc- 

tively with the nonlincarities and increase the order of 
the system In chaotic Hamiltonian systems, investi- 
gations focused on the effect of noise on anomalous trans- 
port and, very recently, on chaotic scattering 0-Q. 

Chaotic scatte ring is a basic process of Hamiltonian 
dynamics B, Hol . with fundamental applications in 
classical [l 2i ] and quantum [HEH systems, and recent ap- 
plications ranging from plankton populations Q to blood 
flows [3] and even the origin of life [HI . Chemical (dis- 
sociative) reactions of simple molecules are also scatter- 
ing processes where chaos is essential in the microcanon- 
ical phase-space formulation of transition state theory 
(TST) jlq |. Scattered trajectories perform transiently 
chaotic motion while trapped by fine structures of the 
phase space, such as fractal nonattracting sets and chains 
of Kolmogorov- Arnold-Moser (KAM) islands @, [ul. 
Noise destroys the small scales of these structures [7J, 
modifies the temporal decay of trajectories in time from 
algebraic to exponential with an exponent that in- 
creases with noise [f| , and create otherwise forbidden es- 
cape paths Q . All these effects weaken the deterministic 
trapping. 

In this Letter we show that noise also plays a construc- 
tive role in chaotic scattering, enhancing the trapping of 
trajectories. First we introduce and scrutinize two dif- 
ferent mechanisms responsible for this surprising effect, 
arguing that they exist in very general circumstances. 
The first mechanism acts in fully chaotic systems and 
reduces the escape rate of particles by blurring the nat- 
ural measure of the system. The second mechanism acts 
on mixed-phase-space systems and enhances trapping by 
throwing trajectories inside KAM islands. We confirm 



the generality of these mechanisms through simulations 
in the conservative Henon map, and we explore the im- 
plications of our results to physical systems. 

To illustrate how noise enhances trapping in fully- 
chaotic systems, consider the baker map [2j 



M 



(x t+ i,y t+ i) = (x t /2,2y t ), 



y < 1/2, 



(x t+1 ,y t +i) = {{x t + l)/2,2y t - 1), y > 1/2, 



defined in [0, 1] x_[0, lb Escape is introduced through an 
arbitrary leak I [l7H2(|: 



M(x t ,y t ) 



M(x t ,y t ) \i{x u y t )£I 
escape if (xt,yt) € /. 



(1) 



For concrctcncss, consider / to be a vertical stripe (/ = 
[x c — A x ,x c + A x ] x [0,1]) at the center of the map 
(x c = 0.5, A x = 0.05). The survival probability P(t) of 
typical initial conditions decays asymptotically as P{t) ~ 
g-*/ 1 " U[l(J 11 1 . In chemical reactions P(t) corresponds 
to the reactant lifetime distribution, which plays a cen- 
tral role in TST [l6| . The trapping strength is quantified 
through the characteristic life time r, which is the recip- 
rocal of the escape rate and is different from the mean 
escape time [3]. It can be obtained from periodic orbits, 
e.g. by calculating the leading root X s of the (truncated) 
polynomial approximation of the zeta function [lfj| , 



1/cm = n( i 



7A P ), 



lnzt, (2) 



where the product is taken over all periodic orbits p, 
which have period t p and expan ding eigenvalue A p , and 
£ = indicates absence of noise |2l|. We locate the or- 
bits analytically so that Eq. ([2]) yields t^ = o with a higher 
precision than simulations. For the baker map exam- 
ple above, all 1990 orbits up to period 18 yield t^ = o = 
6.06 ±0.02. 
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(a) £=0 



(b) 0<^ m <l ( C ) 5=1 



FIG. 1: (Color online) Noise-enhanced trapping for fully 
chaotic systems, (a) The life time r for the baker map |l} is 
shown as a function of the noise intensity £ for two BC: peri- 
odic (red ■) and open (blue •). Horizontal lines: T{=o = 6.06 
(solid) and t* = 9.49 (dashed). Inset: P(t) for £ = 0,£ = 
£ m = 0.035 (open BC) and £ = 1 (periodic BC). (b) Val- 
ues of t*,T£ = o and r m for different positions x c of the leak 
(A x = 0.05). Noise-enhanced trapping occurs when r* > T£ = o 
for periodic BC and when r m > T{=o for open BC (blue shad- 
ing)- 



Next wc consider the effect of noise, added indepen- 
dently to each trajectory x > x+£5 x , y > y + £,5 yi where 
£ controls the noise intensity and 5 x , y G [— 1, 1] are inde- 
pendent uniformly identically distributed stochastic vari- 
ables. Let us first consider the simplest case of periodic 
boundary conditions (BC) for trajectories driven by noise 
outside [0, 1] x [0,1]. This would be the natural choice 
for maps defined on the torus. In this case, for £ — > oo 
trajectories are uniformly distributed in [0, 1] x [0, 1] and 

T^T* With EMI 



= [-ln(l- M (/))]- 



!///(!) for small p(I)}, (3) 



where fi(I) is the Liouville measure (phase-space area) 
of the leak /. This is equivalent to the statistical mi- 
crocanonical estimation of chemical reaction rates in 
TST [jq . In the baker's map example above p{I) = 
2A X = 0.1 and t* = 9.49 which is greater than the noise- 
less case T£ =0 = 6.06 calculated above. The trapping is 
enhanced by noise! 

Let us take a closer look at the two crucial points that 
lead to this simple yet surprising result. The first crucial 
point is the periodic BC that guarantees that r — >• t* 
for £ —> oo. Another natural choice is open BC, in which 
case trajectories outside [0, 1] x [0, 1] escape. This type 
of escape prevails for large £ and reduces r [for the baker 
map r = -1/ ln((l - A x )/(2£) 2 ) ~ l/ln£ for £ > 1], but 
it is negligible for small £. The results shown in Fig.QJa) 
confirm the existence of noise-enhanced trapping for open 
BC, with a maximum r m > T£=o at £ = £ m > 0. Noise 
increases r by more than 50% in the periodic BC case 
and by almost 10% in the open BC case. 

The second crucial point in the derivation above 
is t* > T£ = o- References [17 H20| show that this holds 
for most leaks, also for map ([T]) [18[, what is confirmed 
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FIG. 2: (Color online) Noise tends to uniformly distribute 
surviving trajectories. Quasi-invariant density pj of map JT]) 
obtained at t — 40 for: (a) £ = (r = T£ = q), P5=o distributed 
along the unstable manifold of the chaotic saddle [3. HoL [TH] : 
(b) £ = £ m (r = T m , open BC), pj is smoothed in I and 
negligible outside [0,1] x [0,1]; (c) £ = 1 (t = T*, periodic 
BC), P5=i = Pn is uniformly distributed. The leak / (x c = 
0.5, A = 0.05) and its forward iteration M(I) axe indicated. 



in Fig. [lib). For periodic BC the condition t* > r^ =0 
is sufficient for the existence of noise-enhanced trapping. 
For open BC the condition t* > t^ = o is necessary but 
not sufficient, and noise-enhanced trapping exists only 
when r m > T£ = o- I n Fig. [ljb) this is seen only in the 
deep "valleys" of the t^ = o(x c ) landscape. But the noise 
enhanced mechanism is also present elsewhere, e.g. for 
several x c 's the t(£) curve has local maxima. The dif- 
ference between r* and t^ = q is a consequence of the dif- 
ference between the Liouville invariant density p^ (uni- 
form in x, y) and the open system's quasi-invariant den- 
sity /?{ = o (nonuniform in x, y) inside / [13, El- Note that 
escape occurs one iteration after trajectories enter / (or 
leave [0, 1] x [0, 1]). When p^ = o is large in /, the open 
system measure of / is larger than fj,(I) and t* > t^ = q, 
see Eq. ([3]). In this case the dominant effect of noise is 
to move trajectories outside /, avoiding their escape and 
increasing r. More generally, noise acts on x,y coordi- 
nates making p^ more uniform in x, y @, 0], departing 
from p£=o and approaching p^. This r-increasing effect 
is shown in Fig. |2"1 [TRI^, 

The mechanism described above acts on transient 
chaos and is different from stochastic and coherence res- 
onance|4j . It has been observed in dissipative maps nea r 
crisis [2 




1301 ] and on extended excitable systems [31 
Strong nonlincarity was shown to be a sufficient condi- 
tion for its occurrence in 1-d maps [30j j . No effect was 
observed in a 2D system [32j], raising doubts about its 
generality in higher dimensions. Here we show the exis- 
tence of this mechanism in (Hamiltonian) scattering sys- 
tems. We identify r* > t^ = q as the crucial condition, 
clarifying the generality also for higher dimensions. A 
similar approach in ID maps related this mechanism to 
the nonuniformity of the invariant density of the closed 
system (3^]. In Hamiltonian systems the closed system 
density p^ is uniform, but we show that the mechanism 
is nevertheless effective because the open system den- 
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sity P£=o is non-uniform. 
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FIG. 3: (Color Online) Random-walk model for trajectories 
inside KAM islands, (a) Illustrative KAM island [Henon 
map Q with k — 2], Inset: Magnification around the pe- 
riod three orbit A showing the intersection of manifolds, (b) 
Random-walk model with probabilities p± = 0.5 of stepping 
Ar = ±1, one reflecting (r = 0, center of the island) and one 
absorbing boundary (r = r* ~ r* /£, where f * is proportional 
to island width and £ to step length), (c) P(t) for a simulation 
of 10° walkers started in r = r* — 1, with r* = oo (thin red 
line) and r* = 100 (thick black line). Dashed line: Scaling 
a = 0.5. Inset: Linear-log plot emphasizing the exponential 
tail (for r* = 100) starting at t e « I0 4 = (r*) 2 ~ (r*/£) 2 . 

In mixed-phase-space systems the exponential decay 
discussed so far is replaced by an algebraic decay P(t) ~ 
t~ a due to the stickiness of KAM islands 0,53]. Figure[3] 
illustrates how noise enhances trapping in this case. If £ is 
larger than the lobes of manifold intersections close to the 
KAM island (see inset), it will be more likely for trajec- 
tories to approach (or enter) the region corresponding to 
the KAM island in the deterministic dynamics by jump- 
ing over the manifolds that shield the island. The time 
scales of this effect can be estimated following Ref. @, 
wich argues that for small £ the deterministic algebraic 
decay will be interrupted at a time i& ~ l/C' 3 , where 
P ps 1 depends on a. An asymptotic exponential decay 
after such cut-off is the typical effect of noise in intermit- 
tent systems observed also in scattering systems Q. 
Our approach differs from Refs. 0, HJ because we con- 
sider an additional trapping regime dominated by trajec- 
tories that jump inside the island, circle the elliptic fixed 
point, and perform a random walk in the perpendicular 
direction. This ID random-walk model can be solved an- 
alytically in infinite domain |27[ and yields P(t) ~ t~ a 
with a = 0.5. This is smaller than the deterministic 



exponent 1 < a^ = o < 2, meaning that trapping is en- 
hanced. In reality, the island has a finite area that cor- 
responds to a reflecting boundary in the model, as ex- 
plained in Fig. [3fb). This introduces a cutoff t e > tb 
in the a = 0.5 algebraic decay followed by an exponen- 
tial decay, as shown in Fig. EJc). This second cutoff oc- 
curs when trajectories explore the full island of size r*, 
and since this is a diffusive process r* = r* /£ ~ y/t and 
t e ~ l/^ 2 Q. For £ — > the interval of enhanced trap- 
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FIG. 4: (Color online) Noise-enhanced trapping in the fully 
chaotic (k = 6) Henon map Q. (a) Phase space with trapped 
region, delimited by the solid black line (first intersection of 
the manifolds W u ' s of the fixed point ■). Escape (Inflow) 
denotes an unbounded forward (backward) invariant set that 
diverges in forward (backward) iterations [25|. Initial condi- 
tions and the region where trajectories are removed were cho- 
sen inside these sets [34| . The exit set I is the set of points that 
exit the trapped region in one iteration [24| . (b) Dependence 
of t on £ with T£ =0 = 1.557 (dashed line) and r* = 3.239 
(dotted line). Inset: P(t) for £ = 0, £ = £ m , and with t*. 



ping A< = (t e — tb) —> oo because of the different scalings 
of tt and t e . Therefore, even if not valid asymptotically, 
enhanced trapping is in practice dominant for small £. 
Similar models have been used to investigate different 
problems: anomalous transport Q and escape of initial 
conditions inside KAM islands in random maps [§]. 

Finally, to confirm the general validity of the two mech- 
anisms discussed above, we investigate the paradigmatic 
Henon map (2342^ 



x t+ i = k 



X t -1 



(4) 



where different random perturbations 8t, defined as be- 
fore, are applied to each trajectory. To investigate the 
mechanism for fully chaotic systems we choose k = 6 
in Eq. (|3]), which is the smallest integer for which a 
complete horseshoe exists [26|. In Fig. UJa) we iden- 
tify in the phase space of map (0| the trapped re- 
gion and exit set / (reactant region and transition state 
in TST 0). We estimate graphically fi(I) = 
area(/)/arca(trapped region) = 0.2656, which leads 
through Eq. (3} to t* = 3.239. We compute ana- 
lytically [26j all 226 periodic orbits up to period 10 
and through Eq. © we obtain r ?=0 = 1-557 ± 0.001. 
Again t* > r^ = o, fulfilling the necessary condition iden- 
tified above for the existence of noise-enhanced trapping. 
Our numerical simulations reported in Fig. |U[b) confirm 
the nonmonotonic dependence of r on £ with r m > t^ = q, 
in perfect analogy with the results for the baker map 
shown in Fig. [1] To investigate the mechanism for mixed 
phase-space systems we choose k = 2 in Eq. The 
phase space is similar to the one shown in Fig. U£a) ex- 
cept for the KAM island in Fig. [3)Ja) around the fixed 
point (black dot #). The results are summarized in 
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FIG. 5: (Color online) Noise-enhanced trapping in the Henon 
map with mixed phase space (k = 2). (a) P(t) for £ = 
(thick black line) and £ = 10~ 3 (thin blue line) [H. The 
times tb and t e indicate the beginning and end of the mech- 
anism described in Fig. [3] Power-law scalings are shown as 
a reference (dashed lines), (b) Dependence of tb and t e on £ 
with fits of the model predictions (with /3 = 0.92; see text). 



Fig. [5] and confirm quantitatively the predictions (scal- 
ings of tb, e and a = 0.5) of the random-walk model in 

Fig. EM. 

In summary, we have shown that weak noise leads to 
a slower decay of the survival probability P(t) in fully 
chaotic and in mixed-phase-space scattering systems. In 
both cases P(t) changes nonmonotonically with noise in- 
tensity £ (Figs. @]and[3]), an effect previously observed 
in the dependence of the diffusion coefficient on £ 0, [35| . 
Our approach to the mechanism for fully chaotic sys- 
tems extends naturally to higher dimensions, while the 
mixed-phase-space mechanism has to be expanded to 
take Arnold's diffusion into account 

All scattering systems are likely to be subject to small 
noiselike perturbations and to experience an enhance- 
ment of the trapping. Experiments on the chaotic ad- 
vection of passive tracers [12[ can provide a direct test of 
our predictions: by systematically changing the proper- 
ties of the tracer and fluid one can control the molecular 
diffusion (~ £ 2 ) and obtain the dependence of P(t) on £. 

Our results also describe the effect of noise in the 
phase-space formulation of TST, where an increase of 
the chemical reaction rate is expected. Simulations of 
unimolecular reactions with mixed phase space found 
P(t) t~ a with a < 1 and an exponential tail. According 
to Ref. [l6| the dynamical origins of this behavior remain 
obscure. We show that noise leads to the same observa- 
tions. It remains to be shown whether this explains the 
previous simulations, where noiselike perturbations could 
originate from round-off errors or from higher degrees of 
freedom (see Ref. Q). Our results also impact quantum 
systems. For instance, the increase in the characteristic 
life time r indicates a modification of the fractal dimen- 
sions of the invariant sets [H 0, S| and consequently of 
the fractal Weyl's law 0. 
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